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ABSTRACT 


A long-period forecast is made utilizing a two-level quasi- 
geostrophic model. The model includes friction and heating which is 
a linear function of y. The model was further simplified by restrict- 
ing the disturbance to one wave in the zonal direction. Experiments 
were performed with two dueeatees between the walls. In the case of 
the longer separation, a solution with appreciable time fluctuations 
was obtained with the largest fluctuation being investigated in more 


detail. The smaller separation revealed a traveling wave of constant 


amplitude similar to that observed in certain dishpan experiments. 
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I. INTRODUCTION 


Knowledge of the general circulation has been made possible with 
the advent of high-speed large-capacity electronic digital computers. 
In a classic experiment, Phillips [1956] performed a long-period 
numerical forecast with a two-level quasi-geostrophic model in a 
zonally periodic channel on a 8 plane. Friction and heating were 
included in his equations with the latter as a linear function of 
latitude. He introduced a finite amplitude initial disturbance into 
a baroclinically unstable zonal current. Major deficiencies were the 
error from the 8 plane and the geostrophic approximations, the latter 
making it impossible to account for non-geostrophic dynamics. Further- 
more, due to computational errors, he did not achieve a statistically 
steady state. However, Phillips did initially demonstrate the feasi- 
bility of numerical simulation of the atmosphere. 

General circulation experiments have progressed in sophistication 
and complexity since Phillips’. Smagorinsky [1963] formulated a two- 
level baroclinic primitive equation model by using the Eulerian 
equations of motion within a spherically zonal strip. Continued 
sophistication of the experiments resulted in a nine-level model 
utilizing the primitive equations of motion [Smagorinsky, Manabe 
and Holloway, 1965]. This study included a resolution of boundary 
layer fluxes and radiative transfer involving ozone, carbon dioxide 
and water vapor. A continuation of this study was performed for the 
tropics which was successful in simulating the tropical convergence 


zone and the thermal structure [Manabe and Smagorinsky, 1967]. Mintz 


[1964] conducted an experiment in which the primitive equations of motion 
were developed for a two-level model of the atmosphere with the heating 
and friction terms retained. This model further brought in the effects 
of land-sea contrasts in heating. Nitta [1967] utilized a twenty-level 
model for computing vertical distribution of the geopotential flux and 
conversion of eddy potential energy and eddy kinetic energy due to non- 
stationary disturbances. A model which has the distinguishing feature 

of using height rather than pressure was developed [Kasahara and 
Washington, 1967] and hydrostatic equilibrium was maintained in the 
system. Some advantages of this model are the prognostic equations 

have a simpler form than those of the o system and the lower boundary 
conditions may be easily formulated. The integration of the nine-level 
primitive equation model was further extended by the box method [Kurihara 
and Holloway, 1967] which conserved the conservational properties of the 
original equations. Recently, the effect of hydrology of the earth's 
surface was incorporated in the nine-level primitive equation model 
[Manabe, 1969]. The scheme involved the prediction of water vapor in 

the atmosphere and the prediction of soil moisture and snow cover. The 
numerical integrations were performed for the annual mean distribution 

of solar insolation. 

The increased sophistication and complexity of recent models have 
also produced complex results which make interpretation difficult. 
Therefore, in this study, the complexities have been kept to a minimum 
while still maintaining a description of the basic dynamic processes. 
When the effects of heating and friction are included, the two-level 
quasi-geostrophic model is the simplest model capable of describing the 


general circulation. 


The model is further simplified in that the disturbance is restricted 
to a single wave in the zonal direction. Phillips’ [1956] experiments 
and Smagorinsky's [1963] experiment which contained no forcing function 
in the east-west direction, displayed a single predominate wave number 
which contained much of the disturbance energy. The y-structure of the 
disturbance and the mean flow are calculated accurately with a sufficient 
number of grid points. Some general circulation experiments have been 
made with even simpler models which represent the y-variation of all 
quantities with a few terms of a Fourier series [Bryan, 1959; Lorenz, 
1962, 1963; Young, 1966]. These models do not provide an adequate 
description of the barcerspic interaction between the disturbance and 
the mean flow. Consequently, this model appears to be the simplest 


model which can describe the main mechanisms of the general circulation 


of the atmosphere. 


Ii. THE FORECAST EQUATIONS 


A simple two-level model is constructed by dividing the entire 
ial. AP 
atmosphere into four layers of constant pressure differential, 9 
(fig. 1), numbered 0 to 4 from top to bottom. Assume vertical motion 


to be zero at the top of the atmosphere, while 


0 a 0 

: 1 

Ap ae 2 
. 3 

2 ——— 


Fig. 1. Two-level model used for prediction 


= 


the vertical motion term at the earth's surface is approximated by 


ed 


w= -pEBWw, . (225) 
Charney and Eliassen [1949] used the Ekman theory to derive an 


expression for Wy which is 

1 Z A, Me 

"a 1a Se Aare (2.2) 
O 
where AL is the vertical eddy viscosity and qa, the inflow angle. The 
surface geostrophic vorticity is approximated by 
~ 2 

CA = C3 =V Ys - (2235) 

Following the development of Thompson [1961], but with different 


notation, begin with the quasi-geostrophic vorticity equation 


5 2 2 oY cry 4 
Be YY + kee WYVOV YB, 5 > fo Sea, vem 
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Apply this equation at levels 1 and 3 giving 


j oy Ud) 
Oo vee 2 1 Zz 4 
Be ¥, tik x vy, -VvV WD) Sete: eS ip An Ye (255) 
oh (w, -w, ) 
fom) =2 | 2 3 (hi aaa 4 
me {5 + ik x VY, -V(V ¥3) + a ee an I = AY V ¥,- (2.6) 


If it is desired to stop at some intermediate point, say the tropopause, 


smaller layers could be utilized rather than toinclude the entire atmosphere. 


Next consider the quasi-geostrophic first law of thermodynamics in 





the form 
Qo. oY (S) om _ §a*_ 
ye 5 + lk x VY-V = 1 a ~ Ep O: (2.7) 
where 
2 
RT. or. % 
o= z = C A (2.7a) 


and T. is the temperature from the standard atmosphere. When this 
equation is applied at level 2, the result is 


CA) Ap ow 
veo) i. 
2 alee ates we Gee OO 359) 


oO O 


= (¥)+¥,) +tihkxV 


Define the following quantities 


Yee 

aii 3 
ui aa 9) 3 (2779) 

| 

: eee 3 
Yer = (2. 109) 

which implies 

Yy = ua Ss Yep (ail) 
¥3 = Se = Yen (2.12) 


Here Yep is proportional to the layer thickness and is therefore a 
measure of the mean temperature. Using these definitions, add (2.5) 


and (2.6) and divide the result by 2, obtaining 


il 


3 a a 
Shy tik x VY vvy tk x Y, VV) +8. = - Se 


ran vty Quiey 


For the second forecast equation, subtract (2.6) from (2.5) and 
eliminate Wy using (2.8) 

3 De 2 2 2 
xe (V “U) Vn tik x V¥ VV Mi) ¥p tk x V¥L-VW ey 





oY Wy, 5 1Q, i (2.14) 
+ Seay SS i tron =-u 26” AL, V Yop 
where 
; 2 - % 
Ww = 5 (2.144) 
Ap o 


These are the prediction equations for the model. At this point 
we restrict ourselves to a disturbance of wave number k and the mean 


flow. The fields may be defined as follows 


y E(y,t) FAG@Zt) cas kx FBG se) Heim kx, (2 1S 


m 


F(y,t) + C(y,t) cos kx -ap(y,t) sin kx, (2.16) 


tp 


where k is the x wave number, A through D are Fourier amplitudes of the 
ieee nance and E and F are the zonal mean fields. Now substitute 
expressions for ue and Yer into (2.13), separate the various terms, 
ae eee ine all terms with wave number 2k or higher. Equating coeffici- 


ents of the cosine terms gives 


2 2 3 
O27 ork NT SE O°B + OF OD a7 E OF 
So Ga - ae) = ELS Se * ee 


ee a+ p)]- 8, mn - k°) (A-C) + 


2 
A, eu 5 oa +2 * al, (2.17) 
Oy Oy 
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For the sine terms the result is 


3 3 
2 OE ae A OF OC, QE OF 
a ea abe mae: a 
2 (OE OF ¢ 2 
k (a+ Fc)l +5, ae-x ae heel 
4 2 
a, [2 - 2x? SB - xt B |. (2.18) 
3 oy 
y 


Repeating the procedures for (2.14), the equations are 





& (xe 2 a2) = E op, FB SE 


—— D a cia wi a - 
OF yr OF 
oO. Te ee  - — B - 8 Dk+ 
Oy Sue a | O 
a 2 
K ok?) (ac) + 
Sy 
2 
Mies Ce exe eat fall (2.19) 
and 
~ 
5 (3 - Dk? - Du) =k [sc (Ke “+ we ) + ab ae - u’) - 
y 


3 3 
JE d°c OF SA, OF OE |+ 


6, [22 - Bene Dae 4p |. (2.20) 
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Equating coefficients of terms independent of x in (2.13) gives 


2 2 
0 OE koa as Ome 20 ac 
GC os =. S|) Se 2 a 
ae Loe a 
Ag Ay (2 ae 
dy 


Similarly, coefficients indepdendent of x in (2.14) results in 


(al 2) Fk [, SD pC, FB yp Ray 
2 


2 dy af ay sue 

2 OB OD OA 
alec pec ee, D sage 
aay 2 *S,, aan (2.22 
dy ay. 2 H oy a -22) 


where 


dy | (2.23) 


“The prediction equations are now in the form to be used in the model. 
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Ill. FINITE DIFFERENCE SCHEME AND BOUNDARY CONDITIONS 


ticmr mi ite di ttenernee scheme uSsed@™is illustrated below with a 


sample variable N; 


3 7 Gare ; ee (3.1) 
2 
2 = 2 Car - 2N, + Nev ae 
3 
3 : 3 Guys eee Th a N,) : Gy, - 2N,_, 7 No) (Ges) 


where i is the grid index and H the distance between grid points. 

Centered time differences are used for all quantities except those 
involving friction. The frictional terms are computed at time (t-At). 
In all cases the first step was a forward time step. A special time 
step was introduced every 24 time steps to avoid the separation of 
solutions at even and odd time steps. These special time steps used 
the backwards difference method utilized by Matsuno [1966]. Second 
order equations (2.17), (2.18), (2.19) and (2.20) for time tendencies 
are solved by the exact method of Richtmyer [1957]. A modification of 
the Richtmyer method used by Phillips [1956] was necessary for equation 
(2.22) to satisfy the boundary conditions. Equation (2.21) was solved 
by a direct marching process utilizing the boundary conditions. 

The following boundary conditions were established: 

The quantities A through D = 0 


0 OF Oo Ono _ _ 


ES 


Following Phillips [1956] the mean and disturbance vorticity are 


set equal to zero at the boundary in order to make the boundary walls 


smooth. 
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IV. INITIAL CONDITIONS 


A natural initial state would be one of constant temperature and 
no motion. The heating, in time, would build up a baroclinic current. 
No disturbance could develop until the vertical shear of the Zonal 
current reached a critical value. In order to save computer time, the 
initial state was chosen to consist of a zonal current with a vertical 
shear near the critical value [see Thompson, 1961]. The initial dis- 
turbance has a wave length and a y-structure such that when superimposed 
on a baroclinic current with no horizontal shear it would have maximum 
growth rate. 


The initial conditions for the model are 


A = 200 aia 


B through D = 0 


f HW6 
O 


cos 
W 





£ HW6 
O 


E = —-——- cos 


iM W 
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V. RESULTS 


The main experiments were conducted using the physical parameters 


listed in Table I. 





Fh 
ll 


Q 
il 


1 x 1o™ per second 


squared 


cc a RQ > > > DD 
rt 12, 
Ul ll Ul ll il ll il ll 


x 
ll 


07 oo per meter per second 
290 centibars 

40 meters squared per second 

ox 10° meters squared per second 
22.5 degrees 

2.49 x iowa per meter squared 
200 kilometers 

0.5 hours 


21/4000 per kilometer 


2.3 meters squared per second squared per centibar 


Table I - Numerical Values of the Physical Constants 


The variation of experiments 2 and 3 from the main experiment is 


tabulated in Table 


Experiment 
Experiment 


Experiment 


Table 


W (Km) Q (Kj per ton per sec) 


Es. 

It 8000 4 
Z 8000 8 
3 4000 8 


II - Variation of the Physical 
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lane Gy - = i 


Or (y - = yw 


ae Gy - i 


Parameters 


* 


* 


The forecast period for all three experiments was 300 days. Figure 
2 depicts a plot of the disturbance energy for the 300 day forecast 
period. The disturbance energy referred to throughout the discussion 
is defined in Appendix C (Eqs. C-10 and C-11). A main feature of Fig. 

2 is that the disturbance energy does not reach a steady state although 
it appears to be statistically stable after the first 120 days. 

Figure 3 displays a time section of the mean wind for the period 
120-260 days which is the main time period of interest for this experi- 
ment. The general increase of the disturbance energy begins to level 
off near day 140 which corresponds to the initial appearance of the 
southern jet. This result is similar to observations of Riehl, Yeh 
and LaSeur [1950]. The period between day 190 and day 200 is of con- 
siderable interest. The disturbance energy has its greatest fluctuations 
during this period and it corresponds to the termination of the southern 
jet. 

Figures 4 - 7 represent the amplitudes and phases of ¥ Yep and wo 
in the form Pcos(kx-§) where P is the amplitude and § is the phase. 
Figures 4 and 5 represent day 190, prior to the disturbance peak, and 


Figs. 6 and 7 correspond to day 196, after the peak. The necessary 
oy 
g 





com@ition fom barokropremanstability is that 6 - change sign some- 


oy 
where in the region [Kuo, 1949]. Calculation at day 190 showed that 


this condition was satisfied in the vicinity of the southern jet. It 
is seen in the lower portion of Fig. 4, between y = 1000 km and y = 


2500 km, that the phase tilts of Pes Y_ and w, are opposite to the 


T Zz 


shear of u,- These phase relationships indicate barotropic amplifi- 
cation and therefore a transformation of energy from the mean flow to 


the disturbance. However, the northern portion of Fig. 4 shows the 
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Fig. 4 - Phase Relationship of the Vertical Mean Stream 
Function (Y ), the Stream Function of the Thermal Wind 
(Yr) and " the Vertical Motion (w, = - 22) for day 
190. rs: 
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800 


same tilt as the shear of iam which indicates barotropic damping: this 
implies an energy transformation of disturbance kinetic energy to mean 
kinetic energy. Figure 5 portrays the amplitude of bee and Yip} clearly 
He is larger and both fields have a mid-latitude maximum. A secondary 
peak of amplitude in Yep is beginning to form which is a product of the 
barotropic instability. Figure 6 shows that all parameters (Y Yipe Ws) 


have a phase tilt which indicates barotropic damping over the entire 


range. The amplitude of be and Y,, corresponding to Fig. 6 are displayed 


AE 
in Fig. 7. Of primary note is the distinct double maximum of amplitude 
which is a result of barotropic instability in the southern region. 
Figure 8 depicts the energy transformations over a period which includes 
the maximum disturbance energy. The first item of note is the nearly 
exact correlation between [P-P'] and [P'-K'] which in an energy balance 
would almost completely offset one another thereby leaving a small 
amount of disturbance potential energy change. It should be noted 
that just before the event of maximum disturbance energy, the barotropic 
damping [K'-K] approached a value of zero. This fact permits the dis- 
turbance kinetic energy to increase at a rapid rate due to the large 
value of [P'-K'] and to continue to increase until [K'-K] attains a 
value large enough to terminate growth. 

Figure 9 depicts the disturbance energy over the forecast period 
for experiment 2. As stated earlier (Table II), this experiment has 
twice the heating rate of experiment 1. A steady state was not 
attained but the disturbance energy became statistically steady. The 
main differences from experiment 1 are that the disturbance energy 


fluctuates more rapidly and that its magnitude is much higher. 
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Fig. 7 - Amplitude of 2 and Yep for day 196. 
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Fig. 8 - Energy Transformation for Period which Includes the 
Event of Maximum Disturbance. 
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A depiction of the disturbance energy over the forecast period for 
experiment ae contained in Fig. 10. This experiment differs from 
experiment 1 in that the disturbance between the walls has been cut in 
half. The heating coefficient was doubled eae, its value in experiment 
1. The disturbance energy reached a steady state ee shortly 
after day 120 and stayed relatively constant to the end of the period. 
Thus we have a single wave traveling through the atmosphere with no 
variation in amplitude. This result is similar to certain dishpan 
experiments with an inner core [Hide, 1953; Fultz, et al, 1959]. 

Figure 11 represents the phases of aes Yer ana Wo at day 200 and reveals 
that the Rar ane cece have similar phase variations. The phase tilts 


indicate barotropic damping since there is a single jet in u, (see 


1 
Big. 13). Since Yer lags ¥ there is a northward heat flux and energy 


transformation from the mean potential energy to the disturbance 


potential energy. The correlation mentioned above gives a transformation 


of disturbance potential energy to disturbance kinetic energy since there 
is a tendency for cold air to sink and warm air to rise. Figure 12 
represents a plot of the amplitudes of us and Yer at day 200. The 


amplitude of ie is greater than Y, and both have a single maximum. 


Ak 


Figure 13 depicts the relationship between u,, which could be con- 


ae 


sidered the surface wind, and u, at day 200. This mean wind field is 


1 
Similar to that observed in the atmosphere, although the gradients are 
much less. At the surface, u is westerly in the mid latitudes and 


easterly at high and low latitudes while at the upper levels the flow 


is westerly with a mid-latitude jet. Figure 14 is a plot of the 


Zo 


amplitude of the mean meridional component of the wind (v,) for day 200. 


The amplitude of v. reveals a northerly component at mid latitudes and 


] 
a southerly component at high and low latitudes which implies an indirect 


Ferrel circulation at mid latitudes and a direct Hadley circulation at 


high and low latitudes. 
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Fig. ll - The Phase Relationship of rae te and W» for day 200. 
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Fig. 12 - The Amplitude of - and tor for day 200. 
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Fig. 13 - The Amplitude of the Zonal Component of the mean (u,) 
and the Mean Surface Wind (uy) for day 200. 
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Fig. 14 - The Amplitude of the Mean Meridional Component of 
the Wind (v,) for day 200. 
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VI. CONCLUSIONS 


A two-level quasi-geostrophic model was formulated for simplified 
studies of the general circulation of the atmosphere. The model was 
further simplified by restricting the disturbance to one wave in the 
zonal direction. 

Numerical integration was computationally stable for a period of at 
least 300 days. In all the experiments the same heating function was 
used which was linear in y and independent of time. The main experi- 
ment with a wall separation of 8000 km produced fluctuating solutions 
with irregular time variations. Attention was focused on a sub period 
of the integration where the disturbance featured a rapid growth followed 
by a rapid decay. Apparently a double jet structure in the mean flow 
became barotropically unstable at the beginning of.the period which 
resulted in the F aNBESuEE of the disturbance. By the end of the 
period, a single jet had replaced the double jet structure. Energy 
transformations verify the roll of [K'-K] .in the evolution of this 
disturbance. 

Experiment 2 featured a double heating rate and results were 
similar to experiment 1 with increased fluctuations and higher energy 
values of the disturbance energy. 

Experiment 3 differed from experiment 1 in that the distance 
between walls was reduced to 4000 km and the heating rate was doubled. 
The disturbance energy became constant after approximately 120 days 
and remained so throughout the forecast period. This situation 


represents a wave of constant amplitude propagating through the 
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atmosphere. The behavior resembles the flow observed in certain dish- 

pan experiments which contained a central core [Hide, 1953; Fultz, et al, 
1959]. Apparently the explanation of the difference in behavior between 
experiment 1 and experiment 3 is related to the fact that for W = 8000 kn, 
at least two energy producing modes of disturbance energy can exist in y, 
whereas for W = 4000 km only one energy producing mode can exist. The 
presence of two energy producing modes apparently leads to fluctuations 

in time but the one modal case produces a steady wave. 

Further studies utilizing this simple model could be performed such 
as longer period forecasting; changes in the form of heating; simulating 
dishpan experiments by setting 8 = O and concentrating the heating near 
the boundaries, and a change in the wave number k to name a few. Any 
or all of these studies could further enhance our understanding of the 


- general circulation of the atmosphere. 


of 


APPENDIX A - COMPUTATION OF VERTICAL MOTION 


The computation of vertical motion (5) is necessary for the energy 
equation. To derive the computational equation for vertical motion, 


stare with 


oh oY UQ 
vet _I — | : 
. “a er ese - se SS sy 
where 
l om 
Y =o (A-2) 
m fy y 
1 OY 
Rigi f Se (A-3) 


Substituting equations (A-2), (A-3), (2.15) and (2.16) into (A-1), 
collecting terms, and again neglecting terms with wave number 2k or 


higher results in 


_ 2 oc k (OE OF | 
oy = Teg (p08 a Saige 


aD | k 3 a F E+ ae 


Sin kx /2+4 (A=e) 


Both the disturbance and the average vertical motion are given by this 
equation, which is the equation used in the model for the computation 


of vertical motion. 
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APPENDIX B - COMPUTATION OF AVERAGE NORTH-SOUTH WIND 


The computation of the average north-south wind (v1) is necessary 


for the energy calculations. The computation for Vv, can be derived 


from Fig. B-1l by the following 








ww. ot _ = i 
ae = Q (B-1) (OE 6) a OD 
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Fig. B-l 
In finite difference form equation B-3 becomes 
@, Tels @,) 
= oe , _ why ne os : 
Maer 1), bp (B-4) 


with vy = 0.0. 


APPENDIX C - ENERGY TRANSFORMATION EQUATIONS 


The following definitions are utilized in the energy transformation 


equations of Phillips [1956] 


Mean Potential Energy pe (als : hy dy (C-1) 
Disturbance Potential Energy Pp! Lay, as dy (C-2) 
= 02 oe oe 
Mean Kinetic Energy K = tw Cu, “ Uy ) dy (C-3) 
ay 2 ay 2 Sa 2 ae 
searance Hae G+ G-@] 
Disturbance Kinetic Energy - Le Se = dy 
(C-4) 


Phillips [1956] derived the following disturbance energy trans- 
formation equations which provide an excellent means for examining 
the extent to which the model contains the physical processes known 


to be important in the atmosphere. 
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The symbolic notation of the form [p-p'] for example, signifies a 
transformation of energy from one form - the first in the bracket - to 
the second form. The example [Pp P'] represents the transformation of 
mean potential energy into the disturbance potential energy. 

Utilizing equations (2.9), (2.10), (2.15) and (2.16), the components 


of the disturbance energy become 


W 
Pi = al (c* + D*) ay (C-10) 
and W 
ee af le (07437407 4D") nCOE ec -* E55 en ME coil) 
O 


Similarly, equations (C-5)through (C-9) become 
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